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ABSTRACT 

We present VLA Hi imaging of two blue compact dwarf (BCD) galaxies, 
Haro 2 and Haro 4, and of the spiral galaxy Haro 26, which is projected some 
22' SW of Haro 4. We also show a map of the CO(l-O) distribution of Haro 2 
obtained with the OVRO millimeter array, as well as derive an upper limit for 
C0(2-l) emission from Haro 4 obtained with the CSO. The Hi data of Haro 2 
reveal that the kinematical major axis lies perpendicular to the photometric 
major axis, indicating that the atomic hydrogen rotates about the major axis 
of the galaxy. This confirms earlier indications based on CCD photometry that 
Haro 2 is a dust-lane dE rather than a dlrr. We propose that the present neutral 
and molecular ISM configuration is due to recent gas accretion or a merger. 
The H I distribution and dynamics of Haro 4 and the neighboring spiral Haro 26 
suggest that they are currently undergoing a tidal interaction, reinforcing the 
notion that interactions play an important role in triggering the star formation 
witnessed in Blue Compact Galaxies. 

Subject headings: galaxies: dwarf — galaxies: individual (Haro 2, Haro 4, NGC 3510 
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1. Introduction 

Blue compact dwarf (BCD) galaxies are 
dwarf irregular systems that are currently 
experiencing a burst of star formation. 
Their appearance is dominated by the star- 
burst, which appears compact and blue, 
hence their name. Their spectra show all 
the characteristics of Hll regions, which 
has led many BCDs to be classified as H II 
galaxies (although, in general, not all Hll 
galaxies are necessary dwarf systems). At 
their current star formation rates (SFRs), 
these systems would exhaust their total 
available gas in considerably less than a 
Hubble time; for this reason, it is generally 
accepted that their star formation must be 
episodic. 

Gas-rich dwarf galaxies play an impor- 
tant role in our attempts to understand 
galaxy formation and evolution. BCDs 
show low metallicities which suggests that 
they are chemically unevolved systems. 
They are more metal-deficient than the so- 
lar neighborhood by factors between three 
and twenty (Kunth & Ostlin 2000), with 
a few objects such as IZwl8 (Searle & 
Sargent 1972) and SBS 0335-052 (Melnick, 
Heydari-Malayeri, & Leisy 1992) reaching 
metallicities as low as 0.025 Z Q . 

This said, most BCDs do not seem to 
be young objects in the sense that they 
are undergoing their first episode of star 
formation. Many reveal older stellar pop- 
ulations underlying the bright young stars 
produced in the current bursts (Papaderos 
et al. 1996; Ostlin 2000; Tosi 2001). This 
poses some interesting questions regarding 

1 Present address: Max-Planck-Institut fur ex- 
traterrestrische Physik, Postfach 1312, D-85741 
Garching, Germany 



the production and possible loss of heavy 
elements. Recent X-ray observations with 
the Chandra X-ray Observatory seem to 
indicate that in starbursting dwarf galax- 
ies such as NGC 1569 and NGC 3077 (Mar- 
tin, Kobulnicky, & Heckman 2002; Ott, 
Martin, & Walter 2003), blow-out of ma- 
terial into the halo occurs that will re- 
move the products of nucleosynthesis in 
stars more massive than 8 M Q and dis- 
perse those across a larger volume. Nu- 
merical simulations have suggested that 
in extreme cases blow-away can occur, 
which removes these products from the 
galaxy proper and enriches the intergalac- 
tic medium (Mac Low & Ferrara 1999). 

It is not clear what causes a BCD to 
initiate a burst or episode of star forma- 
tion, and several mechanisms have been 
proposed, ranging from internal instabili- 
ties to external (especially tidal) triggers. 
Reports that some BCDs are isolated sys- 
tems (e.g., van Zee, Skillman, & Salzer 
(1998), who observed four such objects in 
Hi) suggests that interactions may not be 
the only explanation. But other studies 
(e.g., Taylor et al. 1994; Taylor 1997; Wal- 
ter et al. 1997; Pustilnik et al. 2001b, 
and references therein) have shown that 
the presence of close companions produc- 
ing tidal interactions may be a sufficient 
condition to explain BCD bursts. Further- 
more, recent studies (see e.g., Noeske et 
al. 2001) of large samples of star-forming 
dwarf galaxies that look for faint compan- 
ions give additional support to the hypoth- 
esis of interaction-induced star formation 
in BCDs. 

Most studies to date have focused, 
broadly speaking, on the optical charac- 
teristics of BCDs, with several incorpo- 
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rating Hi observations of sufficiently high 
spatial resolution to resolve the objects 
under study. Despite the fact that gas 
has to pass through a molecular phase for 
star formation to commence, only very few 
authors have reported high-resolution (in- 
terferometric) maps of molecular gas in 
BCDs. Detecting CO in dwarf galaxies 
has proved to be quite a challenge, a fact 
that is frequently attributed to the low 
metallicity in these objects (Gondhalekar 
et al. 1998; Taylor, Kobulnicky, & Skill- 
man 1998; Barone et al. 2000; Meier et al. 
2001b). As a result, only a few dwarfs have 
been clearly detected in CO, and still fewer 
objects have been mapped with millimeter 
interferometers (NGC5253: Turner, Beck, 
& Hurt 1997; Meier, Turner, & Beck 2002; 
NGC1569: Taylor et al. 1999; NGC4214: 
Walter et al. 2001; NGC3077: Meier, 
Turner, & Beck 2001a; Walter et al. 2002; 
IC 10: Walter et al. 2004. in preparation). 
With the advent of the FUSE telescope, 
H 2 molecules have been searched for in 
BCDs, but upper limits for the diffuse H 2 
column density are found to be very low, 
with iV(H 2 )< 10 15 cm- 2 (Vidal-Madjar et 
al. 2000). 

To increase the sample of galaxies stud- 
ied in reasonable detail, we have mapped 
the BCDs Haro2 (= Mkn33 = UGC5720) 
and Haro4 (= Mkn36) in Hi. In addition, 
we have obtained an interferometric map 
of Haro2 in the CO(1-0) line. Haro2 is 
catalogued as a Wolf-Rayet galaxy (Beck, 
Turner, & Kovo 2000; Kunth & Joubert 
1985) and has been extensively studied as 
a key object for understanding star for- 
mation activity in low-metallicity environ- 
ments. It is also one of the few BCDs 
that has been detected in CO(1-0), CO(2- 



1) and CO (3-2) single-dish observations 
(Sage et al. 1992; Israel, Tacconi, & Baas 
1995; Barone et al. 2000; Meier et al. 
2001b). 

This paper is organized in the usual way. 
In Section 2, we describe the Hi and CO 
observations. The Hi distributions and 
kinematics of the galaxies are shown in 
Section 3, where we also present the results 
derived from the CO observations. In Sec- 
tion 4 we discuss the atomic and molecu- 
lar gas contents and distributions of Haro 2 
and Haro 4, and the implications for BCD 
evolution. A summary of our results is pre- 
sented in Section 5. 

2. Observations 

EDITOR: PLACE TABLE 1 HERE. 

Table 1 summarizes the optical prop- 
erties of the two galaxies under study 
plus some relevant data on Haro 26 (= 
NGC3510), which lies at a separation of 
~22' southwest of Haro 4. We mapped 
the Hi in Haro 2 and the Haro4/Haro26 
system with the NRAO 2 Very Large Array 
(VLA). Haro 2 was mapped in the CO (1-0) 
rotational transition with the Owens Val- 
ley Radio Observatory (OVRO) millimeter 
array. Single-dish CO (2-1) observations 
of Haro 4 were obtained at the Caltech 
Sub millimeter Observatory (CSO). In this 
section we describe the procedures used to 
obtain and reduce the data. 



2 The National Radio Astronomy Observatory is a 
facility of the National Science Foundation, oper- 
ated under cooperative agreement by Associated 
Universities, Inc. 
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2.1. HI imaging 

The Hi data were obtained with the 
VLA in its C and D configurations, dur- 
ing two runs carried out in November 1994 
and in May 1995 respectively (see Table 2 
for the VLA observing log). Two pointings 
were used for each run, one centered on 
Haro 2 and the other on Haro 4, but includ- 
ing within the primary beam the galaxy 
Haro 26. We observed with the correla- 
tor in mode 4, resulting in two datasets— 
one with 63 channels at 5.2 km s" 1 reso- 
lution, and the other with 31 channels at 
20.8 km s" 1 resolution. Hanning smooth- 
ing was applied on-line during the obser- 
vations. After inspecting both datasets, 
we decided to work only with the lower 
velocity resolution data because of insuf- 
ficient signal-to-noise at high velocity res- 
olution. Both C- and D-array data were 
combined, leading to a spatial resolution of 
about 14"xl5". 

EDITOR: PLACE TABLE 2 HERE. 

Standard VLA calibration and imag- 
ing procedures were applied, using the 
NRAO astronomical image processing sys- 
tem (AIPS, see van Moorsel, Kemball & 
Greisen (1996) for an up-to-date descrip- 
tion). Data cubes were produced with ro- 
bust weighting (Briggs 1995) in order to 
maximize sensitivity while preserving as 
high a spatial resolution as possible. In 
order to define the extent of the Hi emis- 
sion, we initially averaged four line-free 
channels on either side of the band to cre- 
ate a first approximation of the continuum. 
The corresponding image was subtracted 
from the channel maps, resulting in a cube 
of Hi line emission only. We used this 
cube to conduct a more thorough search 



for 21 cm line emission. Following this 
search, a new continuum image was made 
using all channels deemed free from line 
emission. After continuum subtraction, 
the channels containing line emission were 
CLEANed. In order to build moment maps 
from the data cubes we used a masking 
technique: for each CLEANed cube with 
a beam of around 15", we applied a con- 
volution to a beam approximately twice 
as large, i.e., 30". The smoothed cubes 
were then clipped at a level of 2a, and H I 
emission was identified via visual inspec- 
tion on the basis of spatial continuity in 
neighboring channel maps. A conditional 
transfer (or mask) was applied to the orig- 
inal CLEANed cube using the blanked, 
smoothed data. As a last step Hi surface 
brightness maps and velocity fields were 
created. 

EDITOR: PLACE TABLE 3 HERE. 

The final cubes for Haro 2, targeted 
at the center of the observed field, have 
an rms per channel of 0.37 mJy beam" 1 . 
The cubes for Haro 4 and Haro 26 have 
an rms per channel of 0.39 mJy beam" 1 
at the pointing center, which lies close to 
the position of Haro 4. At the position of 
Haro 26, the correction for primary beam 
attenuation elevates the rms by a factor of 
5, approximately. The observations have 
on average an H I mass detection threshold 
for individual clouds of 4.5 x 1O 6 M and 
1.1 x 1O 6 M for Haro 2 and Haro 4 respec- 
tively. These values correspond to surface 
brightness sensitivities of 2.5 x 10 20 cm" 2 , 
where we have assumed that a signal is real 
if it shows up at the 3a level over at least 
two consecutive channels. As a byprod- 
uct of the H I observations, we obtain use- 
ful information about the radio continuum 
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emission coming from the three galaxies as 
well (see below). The Hi observational pa- 
rameters of the combined C- and D-array 
data are listed in Table 3, where columns 2 
and 3 give the beam size in seconds of arc 
and the linear resolution in kpc, respec- 
tively. Column 4 gives the channel width, 
column 5 the rms noise per channel after 
continuum subtraction, and column 6 the 
noise expressed in units of Kelvin bright- 
ness temperature. 

2.2. CO observations 

We observed the CO (1-0) transition in 
Haro2 with the L and H configurations 
of the OVRO millimeter array (Padin et 
al. 1991) during the 1995-96 season. Ta- 
ble 4 summarizes our observations. The 
array includes six 10.4-m diameter an- 
tennas with half-power beamwidths of 
60" at 115 GHz. Two separate correla- 
tors processed continuum and spectral line 
data: an analog correlator recorded con- 
tinuum bandwidths of 1 GHz centered at 
±1.5 GHz relative to the reference local 
oscillator frequency, while a digital corre- 
lator provided 112 frequency channels in 
the upper sideband, Hanning-smoothed to 
4 MHz (10.4 km s -1 ) resolution. 

EDITOR: PLACE TABLE 4 HERE. 

We calibrated the data within the 
OVRO millimeter array database using the 
MMA package (Scoville et al. 1993). Ob- 
servations of Uranus obtained with sim- 
ilar elevation and atmospheric decorrela- 
tion were used to determine the flux den- 
sity of the bright quasar 1150+497 at each 
epoch; we estimate that the resulting un- 
certainty in our flux scale is ~ 10%. Paired 
observations of 1150+497 interleaved with 



observations of Haro 2 every 30-40 minutes 
were then used to remove the effects of in- 
strumental and atmospheric variations on 
phase and amplitude gains. Observations 
of bright quasars were used for passband 
calibration. 

After editing the visibilities in the 
DIFMAP package (Shepherd 1997), our 
final dataset included 24.3 hours of in- 
tegration and 30 distinct baselines. We 
mapped the data using the IMAGR task in 
AIPS, cleaning a single large region down 
to the rms noise in each velocity channel. 
Maps of the 1 GHz continuum channel in 
the line-free sideband showed no evidence 
of emission (with a 3a limit of 2.3 mJy 
at 112 GHz), so subtraction of continuum 
from the spectral line data was unneces- 
sary. Natural weighting produced a syn- 
thesized beam of 3.3" x 2.6" at position 
angle —84.0°, and an rms noise per chan- 
nel of 9.6mJybeam~ 1 . 

We observed Haro 4 in the CO(2-l) 
transition at the Caltech Submillimeter 
Observatory (CSO) on Mauna Kea, Hawaii 
on 2001 April 8. The beam size of the CSO 
at 230 GHz is ~30", and the main beam 
efficiency is determined from observations 
of Mars to be i] m b = 0.72 (Gerin & Phillips 
2000). The 1024-channel 500 MHz AOS 
spectrometer was used for the observa- 
tions. The system temperature was 630 K 
(t"225GHz = 0.2). In total, 1800 s were spent 
on-source, leading to an rms of T\ — 9 mK 
(Tmb = 12.5 mK) in a 6 km s" 1 Hanning- 
smoothed channel. No signal was detected 
in our observations. 
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3. Results 

3.1. HI morphology and kinematics 
of Haro 2 

This is the first extensive interferomet- 
ric Hi work on Haro 2; previous Hi re- 
sults were restricted to single-dish stud- 
ies. Like most star-forming dwarf galax- 
ies, Haro 2 is Hi-rich. Figure 1 presents 
the channel maps. These show clearly re- 
solved structures that are individually ori- 
ented in a southeast-northwest direction, 
and that show a slight displacement per- 
pendicular to this direction as a function 
of velocity. 

The H I parameters we obtain are listed 
in Table 5. Column 1 gives the galaxy 
name; columns 2 and 3 list the veloc- 
ity centroid and FWHM. We conserva- 
tively take half the channel separation 
(10.3 km s _1 ) as the uncertainty for the 
H I central (systemic) velocity. The error in 
the velocity width is estimated to be of or- 
der a quarter of a channel, or m 5 km s _1 . 
The integrated Hi flux (corrected for pri- 
mary beam attenuation) is listed in col- 
umn 4; we estimate an error of about 10% 
for the absolute flux calibration. The radio 
continuum flux determined from the line- 
free channels is given in column 5, and the 
integrated Hi mass in column 6. Column 7 
gives the observed peak H I column density. 

EDITOR: PLACE TABLE 5 HERE. 

Comparison of our Hi detection with 
published single-dish data yields the fol- 
lowing results. Gordon & Gottesman 
(1981) observed Haro 2 with the 300 ft 
Green Bank Telescope, which had a 10.8' 
beam. From examining their published 
spectra, we deduce a peak flux of 34mJy 



and a line width at 50% of the peak of 
At>5o ~ 100 km s" 1 . The authors pub- 
lish a line width at 20% of the peak of 
Af 20 = 208 km s _1 and a total Hi mass of 
5.2xlO 8 M (using a distance of 21Mpc). 
If we adjust their distance to the one 
used throughout this paper (19.5 Mpc), 
the total Hi mass is 4.2 x 10 s M . In 
a more recent study with the Effelsberg 
100-m telescope, which gives a 9.3' beam, 
Huchtmeier, Sage, & Henkel (1995) report 
Av 50 = 110 km s -1 , Av 20 = 170 km s" 1 , 
and a peak flux of 34.0 ± 3.4 mJy. These 
authors point out that their observations 
of Haro 2 may suffer from confusion with 
a close companion located 4.7' to the east 
(see also Nilson 1973). We inspected our 
data cube around the position of the puta- 
tive companion but detected no Hi emis- 
sion. It is therefore justified to use their 
single-dish data to derive an Hi mass 
for Haro 2, which (integrating the pub- 
lished spectra by hand) comes out to be 
3.6 x 1O 8 M . 

The two single-dish results agree with 
each other within their respective uncer- 
tainties. The single-dish linewidths also 
agree with the value we determine from 
our interferometric data, but the latter fall 
short as far as the integrated flux is con- 
cerned: the VLA recovers an Hi mass of 
order 2.3 x 1O 8 M , i.e., only - 55% of 
the single-dish mass (see Table 5). We 
are aware of the fact that CLEANed maps 
do not always yield accurate integrated 
fluxes, especially when the H I distribution 
in the channel maps is extended on scales 
comparable to or larger than ~ 50% of 
the primary beam (see, e.g., Jorsater & 
van Moorsel 1995). In the case of Haro 2, 
whose Hi emission has an extent of or- 
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der 1' — 2', or less than l/15th of the pri- 
mary beam, we discard this clS du possible 
explanation. Instead, we postulate that 
the VLA maps only reveal the "tip of the 
iceberg" and that the galaxy contains ex- 
tended, low surface brightness gas that is 
mostly missed by the interferometer. Fig- 
ure 2 shows the integrated spectrum from 
our VLA data. This profile is in very good 
agreement with that obtained by Hucht- 
meier et al. (1995) in terms of both the 
velocity width and the slight asymmetry 
towards lower velocity side, where a bump 
appears near 1420 km s _1 . We do not con- 
firm the secondary peak seen by Gordon & 
Gottesman (1981) at 1550 km s _1 ; this fea- 
ture was also not detected by Huchtmeier 
et al. (1995) at Effelsberg. In Figure 2 it is 
also clear that we are missing part of the 
extended emission of the galaxy. 

Figure 3 is a map at 15" resolution of 
the H I surface brightness in Haro 2 super- 
posed on an optical DSS image in grey 
scale. It appears elongated in a southeast- 
northwest direction, much like what can 
be seen in the channel maps. The Hi 
disk dimensions are 1.4' x 0.7' (equiva- 
lent to 7.9 kpc x 4.0 kpc at the distance 
of Haro 2), at a level of ~ 5 x 10 20 cm~ 2 . 
These are estimated along the optical ma- 
jor and minor axis of the galaxy, respec- 
tively, as shown by the two lines drawn 
in Figure 3, but were measured from a 
smoothed H I map that is not shown in 
this work. The H I disk orientation roughly 
coincides with the 1.0' x 0.9' (5.7 kpc x 
5.1 kpc) optical disk, whose major axis has 
a position angle of approximately 135°. 
At least three features appear to be unre- 
solved H I clouds with optical counterparts: 
one to the west of Haro 2, (labeled A, at 



« 2 ooo = 10 h 32 m 27 s , 5 20 oo = +54 24 / 00"), 
one to the east (labeled B, at 0:2000 — 
10 A 32 m 37 s , 5 20 oo = +54°24 , 00"), and one 
to the northeast (labeled C, at 0:2000 — 
10 h 32 m 37 a , 5 2 ooo = +54°24'43"). These Hi 
concentrations may either be dynamically 
distinct structures, or represent only the 
local peaks of a smooth extended H I com- 
ponent that the VLA otherwise resolves 
out. In either case, although each of the 
peaks is detected in at least two velocity 
channels, more sensitive observations will 
be needed to confirm their reality. 

Figure 4 presents the Hi velocity field 
of Haro 2 and reveals the galaxy's complex 
kinematics. There is virtually no velocity 
gradient visible along the major axis de- 
fined by the optical light and Hi distribu- 
tion; a cursory inspection of the channel 
maps (Figure 1) already hints at this be- 
havior. A clear gradient is seen at roughly 
right angles to the major axis, with only 
a weak suggestion of rotation about the 
minor axis to the east. The low veloc- 
ity resolution of our observations unfortu- 
nately prevent us from confirming this lat- 
ter feature. In general, the kinematics of 
Haro 2 inferred from its H I emission is con- 
sistent with the analysis of the inner 30" by 
Legrand et al. (1997) using the Ha emis- 
sion line. 

3.2. The molecular gas in Haro 2 

Figure 5 shows the CO(1-0) total inten- 
sity map for Haro 2, integrated over twelve 
channels from 1378 to 1503 km s -1 LSR. 
We detect a clear arclike feature extend- 
ing from southeast to northwest, contain- 
ing several distinct emission peaks. Taken 
together, these cover the same area on the 
sky as that in which vigorous star forma- 
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tion is seen in the optical, as the overlay 
in Figure 6 shows. There appears to be an 
additional concentration of gas ~ 10" to 
the south, and a possible component just 
as far to the northwest (labeled S and NW 
in Figure 5); these clouds are also seen in 
independent data obtained with the IRAM 
Plateau de Bure interferometer by Fritz 
(2000). 

After clipping the data cube at the la 
level, we corrected for the primary beam 
response and measured the galaxy's to- 
tal CO (1-0) flux to be F co = 24.3 ± 
3.2Jykms _1 . Excluding the outlying 
concentrations of gas, we find F co = 
17.6 ± 1.7Jykms -1 . Assuming a Galac- 
tic CO-to-H 2 conversion factor of X co = 
N H JI co = 2 x lO^cm-^Kkms- 1 )- 1 
(Strong et al. 1988; Hunter et al. 1997), 
we estimate the total molecular hydro- 
gen masses following Sanders, Scoville, & 
Soifer (1991): 



M H ^UBX^M S (^)(JL) 2 (- 



.Mpc^ Vjykm 
(1) 

The observed line fluxes then imply 
M Ha = (7.3 ± 0.9) x 10 7 M Q and (5.3 ± 
0.5) x 10 7 M o for all of Haro2 and for its 
central regions, respectively. Multiplying 
by an additional factor of 1.36 to account 
for helium will give the total masses of 
gas associated with the CO-emitting struc- 
tures. Use of a metallicity-dependent X co 
(Wilson 1995; Arimoto, Sofue, & Tsuji- 
moto 1996) would elevate these gas masses 
by an additional factor of 3-4 for Haro 2's 
Z ~ 0.3 Zq. However, more recent stud- 
ies of M 33 (Rosolowsky et al. 2003) and 
metal-poor dwarf galaxies (e.g., Walter et 
al. 2004, in preparation) suggest that Xco 
does not depend on metallicity. Moreover, 
the results of the study of molecular gas 



in Haro 2 by Fritz (2000) validates the use 
of a nearly Galactic conversion factor to 
estimate M Ha in this particular system. 

Since interferometric data are insensi- 
tive to emission which is smoothly dis- 
tributed on scales of order that of the pri- 
mary beam, we wish to estimate the frac- 
tion of the total line flux our map recov- 
ers. We have therefore corrected the entire 
data cube for the primary beam response, 
convolved it to the half-power beamwidths 
of 55" and 22" for the NRAO 12 m and 
IRAM 30 m telescopes, and measured the 
integrated intensities of our convolved data 
at the pointing centers used for published 
single-dish observations. We find we have 
recovered ~ 56% of the the total CO(l- 
0) emission seen by Israel et al. (1995) in 
their 55" beam, and a mean (for four sepa- 
rate pointings) of ~ 61% of the flux seen by 
Sage et al. (1992) and Barone et al. (2000) 
in their 22" beams. We make no attempt 
"00(1-0) fl p correct for this undetected flux in what 
s_ follows. 

A comparison between the CO (greyscale) 
and Hi (contours, as in Figure 3) is shown 
in Figure 7, revealing interesting common 
features. The lower resolution obtained 
for the Hi relative to the CO prevents 
a detailed comparison in the very center 
of Haro 2 where the CO displays a large 
degree of structure. In spite of this lim- 
itation, we confirm that the arc-like CO 
feature follows a distribution along the 
same (southeast-northwest) direction as 
the H I; both are coincident with the major 
axis of the galaxy seen in blue light. The 
CO is more concentrated, with a linear 
size smaller than the stellar component, 
whereas the Hi appears more extended. 
With the exception of the southeast tail, 
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the brightest CO regions are projected 
onto the innermost Hi contour where star 
formation is observed to be most active 
(see also Section 4). 

A more detailed comparison is possi- 
ble between the CO map and the ~ 3" 
resolution RO S AT HRI X-ray image pub- 
lished by Summers, Stevens, & Strickland 
(2001). Figure 8 shows an overlay of CO 
contours on a greyscale representation of 
the X-ray emission. Given that residual 
ROSAT pointing errors (which dominate 
the uncertainty in the relative astrometry 
of these two images) can be as large as 
10", the agreement in the location and spa- 
tial extent of the highest surface brightness 
structures at both wavelengths (as well as 
in the optical; see Figure 6) is quite good. 
In particular, there is a tantalizing align- 
ment between the CO arc and the lower 
of the two X-ray ridges defining a tilted 
"V" morphology in the central part of the 
HRI image. This coincidence can be ex- 
plained if part of the X-ray emission is 
produced by shock-heated gas at the in- 
terface between an expanding superbubble 
and a collection of dense molecular clouds. 
The comparable spatial extent of Ha emis- 
sion in Haro2 (Mendez & Esteban 2000) 
and evidence for an expanding shell in ion- 
ized gas kinematics (Lequeux et al. 1995; 
Legrand et al. 1997) offers further support 
for this picture. Alternatively, the X-ray 
emission could be a consequence of the 
active star formation that is taking place 
currently and the associated detonation of 
core-collapse supernovae. The formation 
of high mass X-ray binaries was deemed 
less likely by Summers et al. (2001), based 
on their hardness ratio as compared to that 
of the extended emission seen in Haro2. 



Higher resolution X-ray observations with 
instruments like XMM-Newton or Chandra 
will be needed to determine conclusively 
the nature of the X-ray emission and its 
relation to the ionized, atomic, and molec- 
ular ISM. 

Figure 9 shows the channel maps for 
the CO(1-0) line in Haro2, obtained after 
smoothing down the data from the original 
10.4 km s -1 resolution by a factor of three. 
Although the signal-to-noise for many of 
the features is rather meager, it is clear 
that the kinematics of the molecular gas 
is quite complex. The central arc-like fea- 
ture in the integrated intensity map (Fig- 
ure 5) does not appear to have a coher- 
ent velocity structure. Moreover, the spa- 
tially outlying structures do not seem to be 
moving at the higher relative velocities ex- 
pected in the case of ordered rotation (e.g., 
the southeast component appears mainly 
in the 1426.6 km s _1 channel). As on the 
larger spatial scales traced by Hi, the gas 
kinematics in the nucleus of Haro 2 are not 
those of a quiescent system. 

3.3. Haro 4 and Haro 26 

Figure 10 shows the Hi channel maps 
for Haro 4. Despite the fact that it is al- 
most twice as close as Haro 2, its Hi dis- 
tribution is only marginally resolved at 
our 14" resolution (equivalent to ~0.5kpc). 
The integrated Hi map superposed on an 
optical DSS image is shown in Figure 11. 
This map confirms the highly compact na- 
ture of the H I distribution, even though it 
is still almost three times more extended 
than the stellar component (see Table 6). 
No clear rotation pattern is detected for 
this galaxy, and higher-resolution VLA ob- 
servations are required to understand the 
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kinematics of this object in more detail. 
The Hi content we derive for this galaxy 
(M ffl = 0.19 x 1O 8 M ) agrees with the 
single-dish value reported by Gordon & 
Gottesman (1981), confirming that no Hi 
emission is missed by the VLA. 

As mentioned in Section 2.2, no CO (2- 
1) emission was detected from Haro4. 
From the per-channel rms of 9mK (T|), 
the corresponding conversion factor 50 Jy K 
for the CSO at 230 GHz (Gerin & Phillips 
2000), and an assumed velocity width 
of 42kms _1 (matching the Hi FWHM), 
we derive a 3cr upper limit on the in- 
tegrated CO(2-l) flux from Haro4 of 

< 21Jykms- 1 . If the CO(2-l)/CO(l- 
0) intensity ratio is ~ 0.5 in temper- 
ature units, as is the case for the two 
BCDs (Haro2 and Mrk297) mapped by 
Sage et al. (1992), then we would expect 
the CO (1-0) line flux from Haro4 to be 

< 43Jykms~ 1 . Equation 1 then implies 
that Haro 4 has a molecular hydrogen mass 
M H2 < 1.6 x 1O 7 M . Correction of X co 
for the possible dependence on metallicity 
according to the prescriptions of Wilson 
(1995) and Arimoto et al. (1996) could in- 
crease this upper limit on Mh 2 by factors 
of 8 and 16, respectively. 

Projected some 22' SW of Haro 4, the 
neighboring spiral galay Haro 26 is well- 
resolved by the VLA beam and spans a 
much larger range in velocity. The chan- 
nel maps are presented in Figure 12. This 
object is a clear example of a galaxy in 
differential rotation seen almost edge-on. 
Hi emission extends from roughly 600 to 
800 km s _1 . Figure 13 shows a compos- 
ite plot of integrated H I profiles of Haro 4 
and Haro 26. The total Hi surface bright- 
ness map presented in Figure 14 shows the 



integrated Hi emission of Haro 26 super- 
posed on a DSS image. As is common, 
the Hi extends well beyond the optical 
image, in this case by about a factor of 
1.5 (see Table 6). We measure a total Hi 
flux for Haro 26 of 40 Jy km s" 1 , in good 
agreement with previous Hi line observa- 
tions. Tifft & Cocke (1988) observed this 
galaxy with the 91-m NRAO telescope and 
obtained an Hi flux of 48.3 Jy km s _1 , 
very close to the 49.4Jykms _1 reported 
by Haynes et al. (1999). More recently, 
Taylor et al. (1994) observed this galaxy 
with the VLA and reported an H I mass of 
10.1 x 1O 8 M (scaled to the our assumed 
distance of 7.9 Mpc), only 15% larger than 
the value of 8.5 x 10 8 M that we report in 
this work (see Table 5). 

Figure 15 shows both Haro 4 and Haro 26 
in a single map at larger scale. Their pro- 
jected separation is ~45kpc. The galaxy 
seen on the DSS image midway between 
Haro 4 and Haro 26 is CGCG 155-052, a 
background galaxy at z ~ 0.03. The 
conspicuous Hi extension to the NE of 
Haro 26, previously reported by Taylor et 
al. (1994), is contiguous with the rest of 
the galaxy in the velocity field (Figure 16). 
The galaxy appears symmetrical only out 
to a radius of ~1' (some 2.3 kpc), beyond 
which the northern tip veers northeast to- 
wards Haro 4. The radio continuum map 
of Haro 26 (Figure 17) shows emission co- 
inciding with this northeast H I feature. To 
the south, isolated radio continuum blobs 
extend farther than the Hi, with intensi- 
ties at the ~ 2.5a level. 
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4. Discussion 
4.1. Haro 2 

Haro 2 is a peculiar member of the BCD 
class. Most BCDs are gas-rich, dwarf ir- 
regular systems manifesting bursts of star 
formation. As first pointed out by Loose & 
Thuan (1986), however, Haro 2 has an R 1 ^ 
rather than exponential surface brightness 
profile in the optical (see also Cairos et 
al. 2001), leading those authors to clas- 
sify Haro 2 as a dwarf elliptical. Haro 2 
has a small H I-to-optical diameter ra- 
tio, Dm/ D opt ~ 1.4, whereas the Hi dis- 
tributions of most BCDs are usually con- 
siderably larger than their optical disks 
(Dm/D opt ~ 2 - 5: Taylor et al. 1994; 
van Zee et al. 1998; van Zee, Salzer, & 
Skillman 2001). The Hi content of Haro 2 
is normal relative to other BCDs, in con- 
trast to its Mm/ L B ratio (see Table 6), 
which is at the low end for BCDs (Smoker 
et al. 2000; Lee et al. 2002; Salzer et al. 
2002). Of the sample studied by van Zee 
et al. (1998, 2001), Haro 2 is most similar 
to UM 38 and Mrk 324, two other dwarf el- 
lipticals (Doublier et al. 1997; Cairos et al. 
2001) with comparable atomic gas proper- 
ties. UM465 resembles Haro 2 in that it 
is a dwarf elliptical detected in both CO 
and Hi emission, although it has a some- 
what smaller H I mass and somewhat larger 
metallicity (Sage et al. 1992). 

While Haro 2 does resemble a BCD in 
terms of its optical morphology and gas 
content, it is almost unique in the kine- 
matics of its neutral atomic and molecu- 
lar gas. Whereas in disk-dominated sys- 
tems the kinematic and isophotal major 
axes are aligned, in the case of Haro 2 they 
are nearly perpendicular. One of the few 



dwarf galaxies with equally pathological 
kinematics is NGC 5253, whose own molec- 
ular and atomic gas components also ap- 
pear to be predominantly rotating about 
the galaxy's major axis (Kobulnicky & 
Skillman 1995; Meier et al. 2002). This 
configuration has been postulated to re- 
sult from a tidal interaction, or from the 
accretion of gas from a formerly gas-rich 
(now stripped) companion. Its kinematic 
complexity, together with the presence 
of a dust lane (like that seen in Haro 2: 
Mollenhoff, Hummel, & Bender 1992), sug- 
gest NGC 5253 may be a dwarf version of a 
giant, post-interaction, dust lane elliptical 
(Oosterloo et al. 2002). In the context of 
this interaction scenario, some of the simi- 
larities and differences between Haro 2 and 
NGC 5253 deserve special attention. 

In terms of star formation, NGC 5253 
appears to be undergoing a burst that is 
markedly younger and more intense than 
that in Haro 2. While the youngest stars 
in Haro 2 have ages ~ lOMyr (Fanelli, 
O'Connell, & Thuan 1988), population 
synthesis fits to star clusters in NGC 5253 
suggest ages in the 1-8 Myr range (Tremonti 
et al. 2001), and its extremely lumi- 
nous and deeply embedded "superneb- 
ula" (Turner, Beck, & Ho 2000; Gorjian, 
Turner, & Beck 2001; Turner et al. 2003) 
must be younger still. Haro 2's older burst 
age is also indicated by its lower millime- 
ter/centimeter continuum flux ratio. In 
particular, our 3a upper limit of 2.3 mJy 
on 112 GHz emission from Haro 2 (see Sec- 
tion 2.2 above) implies that no more than 
3.5 mJy of its 24.6 mJy observed at 1.4 GHz 
can be due to optically thin free-free emis- 
sion (assuming a spectral index a = —0.1). 
Virtually all of the 1.4 GHz flux density 
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from NGC 5253, on the other hand, is ther- 
mal (Turner, Ho, & Beck 1998; Meier et 
al. 2002). Finally, the star formation effi- 
ciency in NGC 5253 is considerably higher, 
with its L m /M U2 ~ 8OOL /M (Turner 
et al. 1997) exceeding the more modest 
~ 5OL /M we infer for Haro2 (see Sec- 
tion 4.3 below). 

It is tempting to associate the older, less 
vigorous starburst in Haro2 with a dy- 
namical disturbance that occurred longer 
in the past than the one currently afflict- 
ing NGC 5253. Several lines of evidence 
indeed favor this interpretation. First, 
much of the molecular gas in the vicin- 
ity of NGC 5253 appears still to be falling 
inward and is not associated with the re- 
gions of most intense star formation (Meier 
et al. 2002). Haro2, in contrast, has a 
CO (1-0) distribution whose spatial extent 
broadly matches that of its star formation 
as traced by Ha, X-ray, and optical con- 
tinuum emission; this agreement suggests 
that Haro 2's molecular component has fin- 
ished coalescing towards the bottom of its 
potential well. Second, the disagreement 
between the position angles of the H I and 
optical disks in NGC 5253- in fact, they 
are nearly perpendicular on the sky- rep- 
resents a level of disturbance not seen in 
Haro 2, where the Hi follows the optical 
disk quite well. Third, we recall from Sec- 
tion 2.1 that our VLA data failed to re- 
cover about half of the Hi emission de- 
tected in single-dish observations of Haro 2 
(see Figure 2). Diffuse, extended emis- 
sion could very well surround the galaxy 
on ~ 10' scales, having accumulated there 
in the same way that late infall in more 
massive post-interaction systems can build 
regular Hi disks (Oosterloo et al. 2002). 



Considering all of the above evidence, 
we conclude that Haro 2 is most likely a 
somewhat older cousin of NGC 5253 — a 
dwarf elliptical galaxy that has recently 
captured gas in a close interaction or 
merger with an initially gas-rich com- 
panion. For Haro 2, an obvious candi- 
date donor is the putative companion lying 
4.7' to the east (Huchtmeier et al. 1995), 
whose non-detection in Hi by those au- 
thors could be explained if it lost all its 
gas in the interaction. If at the same red- 
shift, its projected separation from Haro 2 
would be a mere 27kpc; however, there is 
no redshift estimate for this object, and 
the real distance which separates it from 
Haro 2 remains uncertain. Alternatively, 
Pustilnik et al. (2001b) propose the galaxy 
UGC 5676 as the trigger of the burst of star 
formation. This object lies much further 
away in projection, some 36' (~ 200 kpc) 
northeast of Haro 2, but is known to share 
roughly the same redshift. Both neigh- 
bors are fainter than Haro 2 (by 3 and 1.6 
magnitudes, respectively) and presumably 
less massive as well, consistent with our 
assumption that one of them would have 
lost its gas mass to Haro 2 rather than the 
other way around. In view of the mass- 
metallicity relation for dwarf galaxies (e.g., 
Bender, Burstein, & Faber 1992), Haro 2 
should have received an infusion of less 
enriched gas from its less massive compan- 
ion. This would tend to have left Haro 2 
with a lower metallicity in its interstellar 
medium than in its stars, although exactly 
how much lower is difficult to predict with- 
out knowing what fraction of its gas mass 
originated elsewhere. We note that the 
possibility of tidal interaction with more 
than one object at different times should 
also not be ruled out. 
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In the "passing encounter" scenario de- 
picted above — and even if Haro 2 is actu- 
ally in the final stages of a merger with a 
gas-rich dwarf — gas would have been cap- 
tured at a random angle with respect to the 
major axis of Haro 2, and would have since 
settled into a relatively stable orbit encir- 
cling the galaxy's long axis. The fraction 
of gas settling further into the gravitational 
potential of Haro 2 would accumulate and 
reach densities high enough for molecules 
to form. Molecular gas spiralling inward 
would naturally shock, leading to the cur- 
rent burst of star formation we are witness- 
ing. 

Even if the Hi gas is of relatively re- 
cent acquisition, its motion is by now dom- 
inated by the gravitational potential of 
Haro 2, any plausible donor having passed 
perigalacticon and no longer exerting any 
tidal pull. If we assume that the object 
27kpc to the east was the source of the 
accreted gas, and if we take as an order 
of magnitude estimate a relative velocity 
between the two objects of 150 km s -1 (a 
reasonable upper limit, according to Taylor 
et al. 1994), ~ 180 Myr have elapsed since 
closest approach, i.e., a few revolutions of 
the captured gas in the potential of Haro 2. 
It is therefore valid, as a first approxima- 
tion, to assume for dynamical mass esti- 
mates that the gas is moving in roughly 
circular orbits. We recall that the Hi and 
CO velocity fields show that any rotation is 
within a plane more or less perpendicular 
to the optical major axis. Given that no 
value for the inclination is known, we can 
only estimate a lower limit to the dynam- 
ical mass. The extent of the gas along the 
kinematical major axis is ~1.0'. Taking a 
distance to Haro 2 of 19.5 Mpc and a pro- 



jected rotational velocity of 62.5 km s _1 , 
we find a value of M dyn = 2.6 x 10 9 M Q 
and an Md yn / Lb ratio = 1.0. Even when 
taking into account the fact that the dy- 
namical mass estimate is a lower limit, for 
any reasonable M/L ratio for the stars in 
Haro 2 there doesn't seem to be a need to 
invoke any dark matter. This is contrary 
to what is generally seen in BCDs which, 
if anything, are more Dark Matter domi- 
nated than large, spiral galaxies. 

Our estimate of 180 Myr as the time of 
closest approach leaves ample time for the 
Hi to turn into H2. Following the same 
arguments as given in Braine et al. (2001), 
the time to transform 20% of atomic hy- 
drogen into H 2 is £20% ~ ^ years where 
rim is the Hi volume density (Hollenbach 
et al. 1971, and references therein). A 
rough estimate of the latter can be made by 
taking the average of the H I column den- 
sity across the area where CO has been de- 
tected and assuming a line-of-sight thick- 
ness for the H I of order 200 pc. The H I vol- 
ume density thus derived is rim ~ 3 cm -3 , 
suggesting that the time needed to trans- 
form a substantial amount of atomic to 
molecular gas is less than 10 7 years, much 
shorter than the typical time for solar mass 
stars to form subsequently (~ 30 Myr). 
t 2 o% is an appropriate indicator for the 
timescale on which molecular gas forms in 
Haro 2 as most of the gas is still atomic. 
Lower metallicities and hence lower dust 
content will decrease linearly the rate at 
which this occurs. 

The prominent nuclear star formation 
in Haro 2 has been analyzed in the UV 
and optical by Lequeux et al. (1995) using 
HST. These authors find evidence for gas 
flowing away from the central Hll region 
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at a velocity of 200 km s _1 , notably strong 
Ly a emission with a clear P Cygni profile 
and a complex velocity structure (see also 
Legrand et al. 1997). The column density 
of neutral hydrogen in this outflowing gas 
is N(H) ~ 7 x 10 19 cm" 2 . This is only a 
small fraction of the column density mea- 
sured in the VLA map, which is of order 
2 x 10 21 cm -2 (averaged over the central 
pixel at a spatial resolution of ~ 1.4 kpc). 
The data are best explained by assum- 
ing that the central starburst has created 
an expanding supergiant shell measuring 
~ 1.8 kpc and expanding at 200 km s -1 . In 
this model, the outflowing gas is predom- 
inantly ionized by the central star cluster, 
and the total mass involved in the outflow 
(if distributed in a spherical shell imme- 
diately surrounding the Hll region) is of 
order 1O 7 M . Unfortunately, the size of 
this proposed shell is too close to the reso- 
lution limit of the H I observations for us to 
be able to detect its signature in our data 
cube. 

EDITOR: PLACE TABLE 6 HERE. 

4.2. Haro 4 and Haro 26 

Haro4, shown in Figure 11, appears to 
be highly centrally concentrated in both its 
H I and optical light distributions. The H I 
mass of this galaxy is low compared to that 
of most BCDs previously observed in Hi 
(van Zee et al. 1998, 2001; Pustilnik et al. 
2001b); we estimate a lower limit for its dy- 
namical mass of M dyn = 5 x 10 s M Q . How- 
ever, as no clear indications for rotation are 
seen, this should only be considered indica- 
tive. Other parameters of Haro 4 are rather 
typical among BCDs, such as its M H i / L B 
and Dju/D opt ratios. Kunth et al. (1998) 
find in Haro 4 a broad damped Lya ab- 



sorption profile centered at the wavelength 
corresponding to the redshift of the Ha 
emitting gas. They show that this occurs 
chiefly because the Hi gas is smooth and 
virtually static, as confirmed by the new 
VLA observations, with respect to the ion- 
ized region where Lya photons originate. 

For Haro 26, the Hi map in Figure 14 
confirms the strong warping towards the 
northeast already observed in the optical 
(Sandage & Bedke 1994) and in H I (Taylor 
et al. 1994). While the former authors find 
no nearby companion to explain the warp, 
Taylor et al. (1994) suggest that the emis- 
sion extending towards the north-east may 
in fact be a separate object. We find no ev- 
idence to support this scenario, however, as 
our velocity field in the northeast extension 
appears to be continuous with that of the 
main galaxy. From the channel maps (Fig- 
ure 12) and the zeroth and first moment 
maps (Figures 14 and 16), we speculate in- 
stead that a long-range interaction with 
Haro 4 is responsible for the warp, with 
Haro 4 pulling the near side of the disk of 
Haro 26 towards it. A similar case of a tidal 
interaction between the BCD SBS0335- 
052 and the galaxy NGC 1376 was pro- 
posed by Pustilnik et al. (2001a) to explain 
the burst of star formation in the former. 
The projected separation between the lat- 
ter two galaxies is about 150 kpc, almost 
three times the projected distance between 
Haro 4 and Haro 26. We note that the ra- 
dio continuum emission of Haro 26 shown 
in Figure 17 follows the same warp seen in 
H I. The intensities of the radio continuum 
blobs in the south and the northeast warp 
reach 0.25 mJy beam -1 . 
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4.3. Star formation 

Star formation is thought to occur 
whenever the local gas density exceeds 
a certain threshold. This can be an em- 
pirical level, as proposed by Skillman 
(1987), who found that star formation in 
dwarf irregulars occurs for densities above 
~ 5 x 10 20 cm~ 2 . Alternatively, the local 
gas density can cause clouds to become 
gravitationally unstable (Toomre 1964). 
Taylor et al. (1994) found in their sam- 
ple that this threshold was equivalent to 
a column density of Hi of 10 21 cm -2 . Our 
spatial resolution is insufficient to inde- 
pendently calculate the threshold value for 
Haro2 or Haro4. Qualitatively, star for- 
mation activity is seen in our maps above 
10 21 cm" 2 , and the regions of highest Hi 
column density coincide with the galaxies 
seen in the optical. It is worth noting that 
the CO in Haro2 closely follows the high- 
est Hi concentration, indicating that star 
formation is likely restricted to the very 
central regions as reported by Taylor et al. 
(1994). 

The radio continuum also traces the re- 
cent star formation activity. In the case 
of Haro 2, the radio continuum distribution 
(obtained as a byproduct of the 21 cm spec- 
tral line observations) was only marginally 
resolved and rather featureless. Our 20 cm 
radio continuum flux (see Table 5) com- 
pares reasonably well with the 20 mJy re- 
ported by Beck et al. (2000), who show 
a much higher resolution map than ours 
(VLA A configuration) restricted to a zone 
smaller than the optical galaxy. 

Although Haro 2 is metal-poor, its Hll 
region metallicity at ~ 0.3 Z Q is not ex- 
tremely low. This implies that we are 
not observing its first burst of star for- 



mation and supports the suggestion that 
whereas the youngest burst began less than 
10 Myr ago, it was preceded by other 
bursts (Fanelli et al. 1988). We can esti- 
mate the current rate of formation of mas- 
sive (M > 5M ) stars using the relations 
of Condon (1992). In fact, we can use three 
independent SFR indicators and check if 
they are mutually consistent. 

Based on the 1.4 GHz continuum lumi- 
nosity, and assuming the radio continuum 
is of non-thermal origin (as argued in Sec- 
tion 4.1), we obtain from Condon's equa- 
tion (21) 

SFR1.4 (M > 5M ) _ /. .., 

M yr" 1 ~ 4.0 x lO^WHz- 1 

(2) 

yielding ~ 0.28 M Q yr 1 for Haro 2 if we 
use S'i^ghz from Table 5. Using the Ha 
luminosity, Condon's equation (22) states 

SFR Ha (M > 5 M ) = L Ha 

MqYt- 1 4.4 x 10 34 W [ ' 

yielding ~ 0.39M Q yr~ 1 for Haro 2, if we 
use the total F Ra = 3.8 x 10~ 12 erg cm -2 s" 1 
implied by the narrow-band imaging of 
Mendez & Esteban (2000). Finally, using 
the FIR (40 — 120 fim) dust luminosity as 
obtained by IRAS, Condon's equation (26) 
gives 

SFRfir (M > 5 Mq) _ L FIR 

Moyr- 1 l.lxlO lo L 

(4) 

yielding ~ 0.33 M Q yr 1 for Haro 2. Note 
that we used the IRAS Faint Source Cata- 
log flux densities at 60 and 100 /xm to de- 
rive -Ffir. 3 The good agreement between 

3 We would like to mention that the total SFR in 
Haro 2 must also take into account stars with 
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these three estimates suggests that dust 
does not severely attenuate Haro2's ob- 
served Ha luminosity (at least along our 
line of sight). This favorable geometry, to- 
gether with the disordered CO (1-0) veloc- 
ity field in the inner starbursting region, 
may explain why Haro2 shows a strong 
Ly a emission line with a P Cygni profile as 
reported by Lequeux et al. (1995) and seen 
in new HST/STIS data (Mas-Hesse et al. 
2003). As we emphasized above, the ab- 
sence of similar Ly a emission from Haro 4 
in the same study may be in part related to 
that galaxy's much less disturbed velocity 
field, which renders Ly a photons less able 
to avoid resonant scattering and escape the 
galaxy. 

5. Summary 

VLA Hi observations of Haro 2 reveal 
emission extending over a 7.9 kpc x 4.0 kpc 
area, roughly coinciding with the optical 
galaxy (5.7kpc x 5.1 kpc) in extent and 
orientation. The H I mass we see in our in- 
terferometric data is 2.3 x 10 8 M ; a single 
dish spectrum suggests than that the ob- 
ject is embedded in a smooth, extended H I 
halo containing an additional 1.3 x 10 8 M 
of atomic gas. The M H i / L B ratio is about 
0.1 in solar units. Surprisingly, the kine- 
matical major axis lies almost perpendicu- 
lar to the optical (photometric) major axis. 
This strongly suggest an external source 
for the H I gas. The picture emerging from 
our study is that of Haro 2 being a dwarf 

masses M < 5M Q , the details depending on 
the choice of initial mass function (IMF). Us- 
ing the calibration by Kennicutt (1998) for a 
0.1 - 100 M Q Salpeter (1955) IMF, for example, 
we estimate SFR Ha ^ 1.37 M Q yr _1 and SFR IR ~ 
0.93 MQyr- 1 . 



elliptical galaxy (as already postulated by 
Loose & Thuan (1986) on the basis of their 
optical photometry), which has acquired in 
an interaction or recent merger a substan- 
tial amount of gas. Assuming that this gas 
has settled into regular orbits we derive a 
dynamical mass of 2.5 x 10 9 M & and an 
Mdyn / L B ratio of 1.0. As we have no in- 
dependent handle on the inclination, the 
estimate of the dynamical mass is a lower 
limit and could be larger by a factor of a 
few. Given the uncertainties involved, this 
would imply that there is little need to in- 
voke dark matter in Haro 2. 

High-resolution, interferometric CO(l- 
0) observations with OVRO broadly agree 
with the H I data. The molecular gas coin- 
cides with the region of highest H I column 
density. The orientation of the CO emis- 
sion is aligned with the optical and coin- 
cides with published Ha data. The kine- 
matics, although rather disturbed, resem- 
bles that of the H I velocity field. Assum- 
ing a Galactic conversion factor, we detect 
7.3 x 1O 7 M of molecular hydrogen. 

Haro 4 is unresolved in Hi. We find an 
Hi mass of 1.9 x 1O 7 M . A CO (2-1) ob- 
servation with the CSO yielded an upper 
limit to the mass of molecular hydrogen of 
less than 1.6 x 10 7 M . Haro 4 seems to be 
interacting with the nearby spiral galaxy 
Haro 26 (= NGC 3510), the latter object 
showing a pronounced warp in the direc- 
tion of the BCD. 

Our results provide additional proof 
that interactions with neighboring galaxies 
are responsible for the onset of the current 
episodes of star formation in BCDs. Al- 
though the two galaxies under considera- 
tion lie in relatively low-density environ- 
ments, we conclude that a major accretion 
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event or merger can account for the nuclear 
starburst in Haro2 as well as for its pecu- 
liar kinematics, whereas a more transient 
encounter between Haro4 and the neigh- 
boring spiral Haro26 could explain both 
the star formation in the former and the 
pronounced warp in the latter. 
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Table 1 

Optical properties of the galaxies observed in Hi. 











Morph. 




D 


12 + 




ID 


Alias 


a J2000 


<5j2000 


type 


M B 


(Mpc) 


log[0/H] 


Ref. 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


Haro2 


Mrk33 


10 32 31.9 


+54 24 03 


HII Impec 


-18.1 


19.5 a 


8.4 


1 


Haro4 


Mrk36 


11 04 58.5 


+29 08 22 


BCD 


-14.7 


6.9 6 


7.8 


2 


Haro26 


NGC3510 


11 03 43.6 


+28 53 06 


SB(s)m 


-17.2 


7.9 6 







a Assuming H = 75 km s 1 Mpc 1 . 
b From Tully (1988). 

References. — (1) Mas-Hesse & Kunth (1999) 
(2) Izotov & Thuan (1999) 



Table 2 
vla h i observing log 











VLA 


^int 


Field 


a J2000 


^J2000 


Date 


coring 


(min) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


Haro 2 


10 32 30 


+54 25 00 


1994 Nov 1 


C 


193 


Haro 2 


10 32 30 


+54 25 00 


1995 May 2 


D 


22 


Haro 4/26 


11 05 00 


+29 09 00 


1994 Nov 1 


C 


202 


Haro 4/26 


11 05 00 


+29 09 00 


1995 May 2 


D 


32 
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Table 3 

Parameters of the Hi maps made from the combined C+D array data 





Beam size 


Beam size 


^ ^chan 


rms 


rms in units of Tb 


Field 


(arcsec 2 ) 


(kpc) 


(kHz) 


(mJy beam -1 ) 


(K) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


Haro 2 


14.8 x 14.1 


1.4 


96.6 


0.37 


1.1 


Haro 4/26 


15.5 x 14.0 


0.5 


96.6 


0.39/1.95 


1.1/5.4 



Table 4 

ovro co(l-o) observing log for haro 2 





Array 


Longest 


Passband 


1150+497 


Date 


config 


baseline 


calibrator 


strength (Jy) 


(1) 


(2) 


(3) 


(4) 


(5) 


1995 Oct 5 


L 


115 m 


3C84 


1.0 


1995 Dec 17 


H 


242 m 


0528+134 


0.8 


1996 May 11 


L 


115m 


3C273 


0.9 
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Table 5 

Measured H I and 20-cm radio continuum parameters 



ID 

(1) 


(kms -1 ) 
(2) 


(kms -1 ) 
(3) 


-Phi 
( Jy km s" 

' (4) 


- 1 ) 


5*1.4 GHz 

(mJy) 
(5) 


Mhi 

(1O 8 M ) 

(6) 


» j-peak 

(10 21 cm- 2 ) 
(7) 


Haro 2 


1443 


111 




2.6 


24.6 


2.3 


2.0 


Haro 4 


650 


43 




1.1 


3.4 


0.2 


2.4 


Haro 26 


712 


187 




40 


23.8 


8.5 


2.2 



Table 6 
Inferred global parameters 







Mm/ L B 


M dyn 


Mdyn/ Lb 




M H2 




ID 


Dm/ D opt 


(M /L ) 


(1O 9 M ) 


(M /L ) 


M H i/M dyn 


(1O 7 M ) 


Mm/ Mm 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


Haro 2 


1.4 


0.09 


2.5 


1.0 


0.09 


7.3 


3.1 


Haro 4 


2.8 


0.17 


0.5 


4.8 


0.03 


< 1.6 


> 1-2 


Haro 26 


1.5 


0.76 


11.5 


10.3 


0.07 
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Fig. 1 Mosaic of the channels contain- 
ing line emission from Haro 2, each su- 
perposed on a DSS 5-band greyscale im- 
age. Contours are —2.5, 2.5, 5.0, and 10a; 
negative contours are dashed. The beam 
size indicated in the upper left panel is 
14.8" x 14.1". 

Fig. 2 Global H I profile for Haro 2, ob- 
tained by integrating the channel maps 
over the area of the source after continuum 
subtraction, blanking, and correction for 
primary beam attenuation. For a compar- 
ison with single-dish H I data, we overplot 
the profile (labeled HSH95) obtained with 
the 100-m Effelsberg telescope by Hucht- 
meier et al. (1995) 

Fig. 3 Hi surface density contours for 
Haro 2, superposed on a DSS £?-band 
greyscale image. The contours are 2.0 
(5a), 4.0, 8.1, 12.2, and 16.2 x 10 20 cm -2 . 
The size of the synthesized beam is 14.8" x 
14.1". Three marginally detected Hi 
clouds which seem to coincide with op- 
tical counterparts are identified (A, B, and 
C; see Sect. 3.1). Solid lines indicate the 
directions along which the Hi major and 
minor axes were measured. 

Fig. 4 Intensity-weighted mean velocity 
field of Haro 2. The optical center of the 
galaxy is indicated with a cross; the num- 
bers give heliocentric velocity in kms -1 . 
The synthesized beam is indicated at lower 
left and measures 14.8" x 14.1". 

Fig. 5 Integrated CO (1-0) intensity in 
Haro 2, with offset molecular gas concen- 
trations S and NW indicated. Contours 
are multiples of 0.52 Jy beam -1 kms -1 , for 
the 3.3" x 2.6" synthesized beam shown at 
lower left. Pixels with values of less than 
la were blanked before integrating, and 
the map has not been corrected for primary 



beam attenuation, although the fluxes re- 
ported in the text have been. The cross 
marks the optical position of Haro 2 listed 
in Table 

Fig. 6 Overlay of CO (1-0) contours 
from Figure 5 on a DSS 5-band greyscale 
image. The cross marks the optical posi- 
tion of Haro 2 listed in Table 1. 

Fig. 7 Overlay of H I contours from Fig- 
ure 3 on a CO (1-0) greyscale representa- 
tion (as in Figure 5). The circle at the 
lower left indicates the size of the H I beam. 

Fig. 8 Overlay of CO (1-0) contours 
from Figure 5 on the ROSAT HRI X- 
ray greyscale from Summers et al. (2001). 
The X-ray image has units of count pixel -1 
above background, and has been smoothed 
by convolution with a Gaussian beam with 
a = 3". The cross marks the optical posi- 
tion of Haro 2 listed in Table 1. 

Fig. 9 CO(1-0) channel maps for the 
velocity range integrated to produce Fig- 
ures 5 and 7, after smoothing to 31.2 km s -1 
resolution. Contours are multiples of 
llmJybeam -1 (~ 2a); negative contours 
are dashed. The cross marks the optical 
position of Haro 2 (listed in Table 1); the 
ellipse in the top left panel shows the syn- 
thesized beam which measures 3.3" x 2.6". 
These maps have not been corrected for 
primary beam attenuation. 

Fig. 10 Mosaic of the Hi line emission 
channels of Haro 4 overlaid on a DSS B- 
band image. Contours are -2.5, 2.5, 5.0, 
10, and 15 a; negative contours are dashed. 
The beam size indicated in the upper left 
panel is 15.5" x 14.0". 

Fig. 11 H I surface brightness contours 
for Haro 4, superposed on a DSS £>-band 
greyscale image. The contours are 2.1 
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(5cr), 6.2, 12.3, and 18.5 x HPcm" 2 . The 
beam size is is 15.5" x 14.0". 

Fig. 12 Mosaic of the channels show- 
ing Hi line emission from Haro26, su- 
perposed on a DSS £>-band greyscale im- 
age. Contours are —2.5, 2.5, 5.0, and 10 a; 
negative contours are dashed. The beam 
size, indicated in the upper left panel, is 
15.5" x 14.0". 

Fig. 13 Global H I profile for Haro 4 and 
Haro 26, obtained by integrating the chan- 
nel maps over the areas of the sources after 
continuum subtraction, blanking, and cor- 
rection for primary beam attenuation. 

Fig. 14 H I surface density contours for 
Haro 26, superposed on a DSS S-band 
greyscale image. The contours are 2.1 
(5 a), 6.2, 12.3, and 18.5 x 10 20 cm- 2 . The 
beam measures 15.5" x 14.0". 

Fig. 15 A larger scale Hi map of the 
field including Haro 4 and Haro 26 (at up- 
per left and lower right, respectively), over- 
laid on a DSS B-band greyscale image. 
For clarity, only the two lowest contours 
from Figures 11 and 14 are shown. 

Fig. 16 Intensity-weighted mean veloc- 
ity field of Haro 26. The numbers indicate 
heliocentric velocity in kms -1 . The beam 
size is indicated by the ellipse and mea- 
sures 15.5" x 14.0". 

Fig. 17 The radio continuum emission 
for Haro 26, obtained by averaging the 
line-free channels of the Hi data cube. 
Contours are 2.5, 5.0 and 7.5a (la ~ 
0.1 mJy beam -1 ) and are superposed on a 
DSS B-band greyscale image. The beam is 
indicated in the lower left corner and mea- 
sures 16.8" x 15.6". 
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